o 



Can the Big Bang Singularity be avoided by a single scalar field? 
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In this note, we investigate the possibility of avoiding the Big Bang singularity with a single 
scalar field which couples non-minimally to gravity. We show that in the case that gravity couples 
linearly to the field, some severe conditions on the field's potential have to be imposed. However, 
in non-linear case, it is quite generic to avoid the singularity with single scalar field. 

I. INTRODUCTION 

The Big Bang Singularity (BBS) is a notorious problem that present in Standard Big Bang Theory as well as 
inflation theory \1\. People have made a lot of efforts in solving this problem, and lots of mechanisms have been 
carried out, among which there are scenarios from string theory, quantum gravity theory and so on. Either reduced 
from these theories, or directly arisen from the effective field theory, we are able to obtain a bounce solution in 4- 
dimensional space-time. Bouncing cosmologies, especially non-singular ones, has gained more and more development 
in many aspects |2| ■ As an alternative of the standard theory in the early universe, it can not only avoid BBS, but also 
gives well favored observational effects compared to the astronomical data, such as scale invariant power spectrum 
and suppression of CMB low quardrupole on large scales [3j. Besides, a scenario of bouncing cosmology can give rise 
to large nongaussianities [4|]. 

If we take our universe to be flat, thus within the framework of Einstein's General Relativity in 4 dimensional 

i-G space-time, a bounce is usually realized by Quintom matter ,5] 28]. As has been briefly proved in [5j], in order for a 

'—^ bounce to happen, where the scale factor of the universe contracted to some value and then expanded, the Hubble 

parameter has to evolve from a negative value to a positive one. It will cause the equation of state of the universe will 

D \ be less than —1, violating the Null Energy Condition (NEC). Generally, some Phantom part will get involved that 

i-G • makes the physics unclear Q. 

However, this dilemma can be ameliorated in many ways, such as to introduce non-scalar fields, or higher derivative 

£»0 . operators and so on. One of the simplest ways is to go with the help of non-minimal coupling to Gravity [8j. In 

£> quantum gravity in curved space-time, it is argued that the existence of the non-minimal coupling term is required 

0"N , by the quantum corrections and renormalization [9j. This coupling includes scalar-tensor theories where the Ricci 

Cn i scalar couples to the field th roug h the term F(d>)R UM fas will be called "linear coupling" for convenience), which 

contains Brans-Dicke theory [ll| or dilaton theory [1J, [l3| , see [Tj| for a comprehensive review. Since the inclusion 

of non-minimal coupling term can drive cosmic acceleration with a wider class of potential than usually considered, 



(N 



it can also be utilized as inflaton in the early time [l5| and dark energy at current epoch [lj|. With the help of 



o 
o 



X 



the non-minimal coupling, a single scalar field can also behave like systems of multi-degree-of-freedom without really 
involving a ghost field, such as having its equation of state (EoS) cross the cosmological constant boundary as a 
"Quintom" matter |17| |29l]. M oreover, a great deal of non-minimal coupling theories can find their equivalence to the 
modified gravity theories |20| . 

In this paper we study the possibility of realizing bouncing scenario with a single scalar field non-minimally coupled 
to Gravity. We find that generally speaking, it does work as expected, except that for some specific case where Gravity 
couples linearly, some constraints have to be imposed on the field action. The paper is organized as follows: in section 
II we shortly demonstrate the general case for a non-minimal coupling field to drive a bounce, in section III and IV 
we investigate the linear coupling case and non-linear coupling case respectively, while conclusions come in the last 
section. 

II. GENERAL CONDITIONS FOR A SINGLE FIELD TO AVOID BIG BANG SINGULARITY 

To begin with, let us consider the most general action containing one scalar field and gravity, of which the two 
components are coupled in a very general form: 

s = [d?xV=9{£%fi + P(x,4>)} , (1) 
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where n 2 = 8irG and the metric used here is g^ = diag[l,—a 2 (t),—a 2 (t),—a 2 (t)]. X denotes the kinetic term of 
the field: X = iV^V^. What we're working throughout this note is under Jordan frame [301] . By variation of the 
action we can easily obtain the Friedmann Equations: 



2 

/(£,<£) 



- k 2 P(X, 0) + k 2 P x 2 + 3(# + P 2 ).//j - 3>Hf R = , (2) 

k 2 P(X, 0) ~f R (H + 3P 2 ) + f R + 2Hf R = , (3) 



where /r = t%, and the equation of motion for the field </> is: 

- P X D<£ - P XX V^XV^ - Px^V^V^ + P + A = o . (4) 

where D<^> = + 3P0 and similarly /^ = ^4 is the derivative of /(P, ^) with respect to the field </>. 

As mentioned before, for a flat universe with standard Einstein Gravity in 4-dimensional space-time, a bounce 
happens only when the Null Energy Condition (NEC) is violated. We can learn from Einstein Equations that this 
leads to the vanishing of Hubble parameter H with a positive time derivative, namely H — and H > at the bounce 
point [Sj. Substituting the first condition to the above equations one gets: 

3/r + k 2 (p + 3P) + / = 0, (5) 

H = ^>0. (6) 

In following sections, we will analyze in detail the application of these conditions to the non-minimal coupling scalar 
field. Before starting, we would like to classify the coupling term into two categories according to its properties: one 
is a more specific case that the Gravity couples linearly to the field, and the other is a more general case which refers 
to non-linear coupling. Actually, these two categories are very different due to the different numbers of degree of 
freedom in the background equations. Due to this reason, it may cause very different conclusions. 

III. CATEGORY I: GRAVITY COUPLING LINEARLY TO THE FIELD 

A. general solution 

Among the various possible forms the coupling terms may have, it is the simplest to study the one that Ricci scalar 
couples linearly to the field, which is easy to get a renormalizable theory without any suppression. Moreover, as differ 
from its non-linear counterpart, in this case the higher order derivative of R is zero, thus higher order time derivatives 
of H will not get involved into the background equations. To study the general solution of this case, let us assume 
that /(P, <f>) — P(l — F ((/>)), where the first term in the bracket denotes the standard term from General Relativity, 
and the second comes from the non- minimal coupling. From Eqs. @ and ([3]), one can get: 

R(l-F(4>)) , _ 2 



P(1-P(0)) , , 



+ k' P + 3H(1-F(4>)) =0, (7) 

+ K 2 P(X, 0) - P (1 - P(0))-P(0) = , (8) 



2 
which can be furtherly simplified to a very neat form: 

-3H 2 (l-F) + n 2 p + 3HF = , (9) 

2k 2 XP x + 2P(1 - F(cf>)) - 2XF U - F^ = . (10) 

In deriving these equations, we have made use of the fact p — 2XPx ~ P and X — cfxf). The last one of the above 
equations gives 



2k 2 XP x + 2P(1 - F(<f>)) - 2XF 4> 



(11) 
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For H = 0, we have 



p = , (12) 

* - " 3F* + 2 K > PX {1-F) ■ (13) 

In deriving the last equation we also made use of Eq. (j4]). We can see from above that for case of gravity coupling 
linearly to scalar field, a bounce requires the rhs of the above equation to be larger than as well as the energy density 
equals to 0. This can be viewed as a general condition for this case. 

B. single scalar field with canonical form 

At first let us consider the simplest case in which the field has a standard canonical form: 

1 
2 



P(X,4>) = ~d^<t>-V(cf > ), (14) 



where V{4>) is its potential. Thus the energy density of the field p = ^cp 2 + V(4>). The condition of p = requires 
that at the bounce point \<\? + V(<p) — 0, or equivalently V{4>) = —\<j) 2 - It can be straightforwardly seen that as long 
as 4> 7^ at this point, one needs V{4>) < 0, namely, a potential containing negative values at least at some region 
is required, while a positive-definite potential as usual doesn't work. But is that the case even for the possible case 
where <p = V((f>) — 0, though it refers to a very specific point in (</>, cf>) space and needs strong fine-tuning? To see 
this, let's take a look at the other condition of H > 0, which then indicates: 

where we used p^ = V$ as well as px — 1. From this equation we notice that V^ can be either larger or less than 
0, depending on the form of F(<f>). For example, if we assume that F(<f>) = k 2 £0 2 , then the condition is ^,<j)V^ < 0. 
But anyway, as long as V$ ^ at the bouncing point where V((f>) = 0, there will inevitably be some region where 
V((f) < 0. So this will be a generic constraint on the potential of the field. In usual case, this will break down the 
positive-definition of the energy density, which causes unphysics. However, here it is fine both because of the positive 
kinetic term and the effects of gravity. Thus, after some proper choice of parameters, we can still get a self-consistent 
system. 

As a further study, we also derived the second derivative of H with respect to t at the bouncing point. It appears 
to be: 

H = [2k2(1 _l )+ 3f2]2 {[2k 2 (1 -F) + 3F 2 }Ik 2 (-2 K 2 J>4> - 4ff^ - AH^ - AH^F^ 

-V^F^ - V^F^4> + F^<j> 3 + 2F H ,M>) ~ 12HHF 2 ~ \2H 2 F^F^4>] 

+ (2k 2 F^ - 6F^<J))[k 2 {~k 2 <P 2 - AHQFj, - V4F4, + F^ 2 ) - 6H 2 F 2 }} . (16) 

It can be straightforwardly read off that for the <f> = V((f>) = case (note that H = 0), H will vanish at the bouncing 
point. This means that in this case, the velocity of Hubble parameter always reaches its extreme value when passing 
through the bouncing pivot. This is a new property that hasn't been pointed out by other authors in the literature. 
In order to support our analytical calculation, we also try to find some regions in parameter space which can 
give good numerical results. In order to satisfy the constraints on the potential, we take into account two forms of 
potentials: 1) V((f>) = ^m 2 (j) 2 — Vq and 2) V((f>) = Vo(e x ^ — |). We add some negative zero-point energy to the 
potential to let it get some negative region. For the coupling term, we choose F{4>) = k 2 ^4> 2 , where £ is some coupling 
constant. In Figs. Q] and [2] we show some numerical calculations on Hubble parameter H and scale factor a. We can 
see that a bounce can happen naturally with H crosses the zero divide line. This indicates that by this way, the Big 
Bang Singularity can be avoided. 

C. single scalar field with DBI form 

Beside canonical ones, there is another form of scalar field that is very common used in cosmology: the Dirac-Born- 
Infeld (DBI) form. Actions of this form effectively describes tachyon dynamics, which can be obtained naturally in 
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FIG. 1: The case of canonical single scalar field and linear non-minimal coupling. The potential is chosen to be V(4>) 



Im 2 - 2 



Vb and the coupling term is RF{4>) = «; £-R0 . The parameter is chosen to be m = 0.1m p i, Vb = 0.005m p ; and 



£ = —1.0, and the initial conditions are (f>i = 0.473m p ;, (jn — 0.061m p! . 
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FIG. 2: The potential is chosen to be V(<p) = Vb(e A< * — \) and the coupling term is RF(4>) = n 2 £,R(f> 2 . The parameter is chosen 



to be A : 



1.0m j , Vb = 0.05m ; and ^ = 1.0, and the initial conditions are 4>i 
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string theory [21]. Moreover, it can also drive the acceleration of the universe and can act as inflaton [22j] and dark 
energy [23(. As inflationary models inspired by the string theory, the DBI action non-minimally coupled to gravity 
has also been investigated in the literature, cf. Ref. J24J. Here we study the last case to see whether it can avoid the 
Big-Bang Singularity. 

The DBI lagrangian containing non-minimal coupling term is: 



P(X,< 



■V((j))Vl - 2aX , 



where a is some positive constant. Its energy density can be calculated as p = 2XPx — P = , }t „ 
obtain its equation of state and sound speed square by their definitions: 

- = 2aX - 1 , c 2 s = ^ 
P Px 



w 



Vl - 2aX , 



(17) 
One can also 

(18) 



which shows that this kind of action is stable under classical perturbations (0 < c? s < 1) with Quintessence- like 
behavior (— 1 < w < 0). 
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From (fl7|) we obtain the equation of motion for the scalar field: 



aV 



-□« 



a 2 V 



-Xcj) 



nV " r <p-V*VT=toX-^=0 



where we get: 



(l-2aX)i (l-2aW)i (1 - 2aX)i 

2 _ [2K 2 (3HaVci> + V4,) + RF^l - 2aX)*](l - 2aX) 



2k 2 



2K 2 a^ 



and the Friedmann equation: 

3ff 2 (l-F(0)) = K 2 V ^ ] 



ZHF^ , 



H = 



-k 2 [k 2 2aVX , 

(l-2aX)2 



+ aVHF^tj) + F (l - 2aX)(3.ff0aV r + V*) - 2^X0^] - 6H 2 F 2 {1 - 2aX)i 



2k 2 (1 - F)aV + 3F 2 (1 - 2aX)§ 



(19) 



(20) 



(21) 



(22) 



From above we can see that the first condition at the bounce point p = only requires V = 0, which is looser than 
that of the canonical scalar field case. Moreover, F^Vq < is still required by the second condition H > 0. Therefore, 
in order to have a bounce, a region of negative value of V{<p) is a ls° inevitable. 

We also make the numerical calculations, with the potential of the form V((f>) — Vo(e x ^ — 5) and some proper 
parameter choice. The bounce happens naturally from the view of the plots. As a side remark, one may notice that 
due to different initial conditions, the model may present various behaviors. From Fig. ([3]) we can see that, the hubble 
parameter varies very fast after the bounce, from increasing to decreasing, indicating that the universe will enter into 
a moderate accelerating or decelerating expansion phase soon. While from Fig. Q, one may note that the hubble 
parameter experienced a period of slow variation. This indicates that it is also possible to give rise to an inflationary 
period after bounce, although we will not discuss in detail as it goes beyond current topic. 




FIG. 3: The potential is chosen to be V{4>) — Vo(e — |) and the coupling term is RF(<f>) = n 2 £R(f) 2 . The parameter is chosen 
to be Vq = l.Orripi, a — l.OmL 4 , A = — l.OmZ 2 and £ = 1.0, and the initial conditions are (pi — 0A07m p i, <j>i — 0.776m 2 i. 



IV. CATEGORY II: GRAVITY COUPLES NON-LINEARLY TO THE FIELD 

As a comparison, in this section, we will focus on a more general case where, in the coupling term, Gravity possesses 
a non- linear form. The most general form of basic equations comes from Eq. ((2E|), However, for the sake of simplicity, 
we only discuss about the canonical field (|14[) . Therefore, we can derive the Friedmann equations and the equation 
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FIG. 4: The potential is chosen to be V((p) = Vb(e A< * — |) and the coupling term is RF((f>) = K 2 ^R(f> 2 . The parameter is chosen 



to be Vo = 1.0m pI , a — 1.0m 



"pi 



-1.0m , and £ = 1.0, and the initial conditions are 4>i — 0.500m p ;, 



0.548m^ 



of motion explicitly: 



H 
H 



where 



*, [(Rf R - /)/2 - 3HM + « 2 (^ + V(0)) - 3/^ff 2 ] - 4ffff , 



1 



6/rr 



[ff (/flfli? + /fl00) - 2/flif - K 



2 J.2 



- 3#0 + V 



/</> 



. 



(23) 



/ijflfl^ 2 - 2f RR4> ^R - Jr^4> 2 ~ fR4>4>] - ±HH - 4H 2 , (24) 



(25) 



Since the non-linear term of Gravity has been involved, the higher order derivatives of H appears in the equation, 
which makes the equations difficult to solve analytically. But as the order of derivative increases, the number of 
effective degrees of freedom becomes more, and it will be easier to violate NEC and realize the bouncing process. 

Figs. [5] and |5] show that a canonical field with coupling term contains logarithm function of R. The logarithm 
coupling may seem strange, however, one can find the similar form in previous phenomenological studies, see e.g. [25j. 
We can see from the plot that with proper choice of parameter, this kind of coupling may also give rise to a bounce 
scenario. Fig. [7] is the plot of universe behavior for the coupling of R 2 to some exponential potential of (f>. This kind 
of potential looks like the dilaton potential in string theory [12T |. With this kind of coupling, the universe can also 
pass through the bouncing point smoothly. 



CONCLUSIONS AND DISCUSSIONS 



In this note we investigated the possibilities of a single scalar field giving rise to a bouncing scenario in the very 
early universe, with some non-minimal coupling to gravity. Along with some explicit examples, we showed that it is 
quite possible. However, for the common case that Gravity couples linearly to the field, the field will suffer from some 
severe conditions, such as abandoning the positive-definition of the potential. For other cases where Gravity couples 
non-linearly, it will be more free to get a bounce since more degrees of freedom have been evoked. 

The realization of bounce with non-minimal coupling field is important in the sense that, with the involvement of 
Gravity, one no longer needs any ghost field to violate NEC, which might cause problems. Recently there have been 
many works on construction of bouncing cosmologies within the framework of non-minimal coupling theories, which 
have some relation to our work while focus more on their fundamental origins [26| . Furthermore, it is expected that 
with proper choice of the coupling form, it is also possible to obtain right amount of observational signatures to meet 
the data, and to give rise to some new features and predictions for future experiments, for example, a recent work 
shows that the matter bounce scenario with non-minimal coupling will give rise to scale-invariant spectrum and large 
particle productions [27|. All these fancy topics are left for the forthcoming work. 
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FIG. 5: The potential is chosen as V(4>) — \m 2 <j> 2 and the coupling term is chosen as f(R,cf>) — R — 2K 2 £ln (-^-)(f> 4 . The 
parameter is chosen to be m — 1.0m p i, £ = 1.0 and Rq — 1.0m 2 ;, and the initial conditions are 4>i = 0.783m p j, cj>i — — 0.205m 2 ; . 
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FIG. 6: The potential is chosen as V((f>) — \m 2 § 2 and the coupling term is chosen as f(R,(j>) — R — 2k 2 £ui( 
parameter is chosen to be m = 1.0m p ;, £ = 1.0 and Ro = 1.0m 2 , ; , and the initial conditions are 4>i — 0.500m p ;, <f>i — 0.036m 2 ;. 
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